Validity of numerical trajectories in the synchronization transition of complex systems 
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We investigate the relationship between the loss of synchronization and the onset of shadowing 
breakdown via unstable dimension variability in complex systems. In the neighborhood of the 
critical transition to strongly non-hyperbolic behavior, the system undergoes on-off intermittency 
with respect to the synchronization state. There are potentially severe consequences of these facts 
on the validity of the computer-generated trajectories obtained from dynamical systems whose 
synchronization manifolds share the same non-hyperbolic properties. 



There is a steadfast interest in spatially extended dy- 
namical systems, as coupled map lattices, in which both 
space and time are discrete variables, due to their vari- 
ous physical and biological applications . Complete or 
identical synchronization (CS) occurs when all the cou- 
pled maps have the same value for their dynamical vari- 
ables @. Weakly coupled maps are not usually identi- 
cally synchronized, but as the coupling strength increases 
they can reach CS. The approach to this state occurs in 
an intermittent fashion, alternating quiescent states of 
synchronized behavior with irregular bursts As the 
average duration of these quiescent states tends to infin- 
ity, we attain a CS state. The onset of the transition 
from CS to intermittent synchronization indicates that 
the dynamics of the synchronized state ceases to be hy- 
perbolic through a mechanism called unstable dimension 
variability (UDV) 0. 

A dynamical system whose invariant set (like a chaotic 
attractor) exhibits UDV fails to have the splitting be- 
tween stable and unstable manifolds which is consistent 
with the dynamics. The consequences of UDV are poten- 
tially dangerous from the point of view of the reliability 
of the numerical trajectories one gets from the dynamical 
system 0. A computer-generated chaotic trajectory is 
unavoidably subjected to one-step errors caused by us- 
ing finite-precision arithmetics, which makes the noisy 
trajectory to exponentially diverge from the true chaotic 
orbit it is intended to represent . However, when the dy- 
namical system presents the shadowing property, there 
exists a true chaotic trajectory which stays uniformly 
close to the numerical trajectory for a certain time in- 
terval 0. This shadowing trajectory may not be that 
particular orbit we have been looking for, but for many 
purposes - as when computing statistical quantities - they 
can be equally useful 
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In this communication, we aim to harness the connec- 
tion between the intermittent transition from CS and the 
shadowing breakdown induced by UDV, so as to explain 
the intermittency features from the statistical proper- 
ties of the appropriate finite-time Lyapunov exponents. 
Moreover, if the shadowing breakdown is so severe that 
we could not find a true chaotic trajectory that stays 
close to the numerical trajectories we generate, the lat- 



ter have their validity compromised. As we will see, this 
will be typically the case in the transitional regime from 
CS. Previous papers .8] have dealt with shadowing break- 
down via UDV in lattices of maps and oscillators with lo- 
cal coupling. On the other hand, the present work deals 
with non-locally coupled maps, for which a transition to 
CS is observed as the effective coupling range is varied, 
from a global to a local interaction form 

We examine a one-dimensional chain of N coupled 
chaotic logistic maps at the outer crisis, x <— * f(x) 
t:r: 1 — x), where Xn £ [0, 1] represents the state variable 
for the site z (z = 1,2,... N) at time n, with a coupling 
prescription in which the interaction strength between 
sites decays in a power-law fashion with the lattice dis- 
tance |l(J 
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(1) 

where < e < 1 and a > are the coupling strength 
and effective range, respectively, and r](a) = 2 XwLi 3~ a > 
with N' = (N - l)/2, for N odd. This coupling pre- 
scription allows one to pass continuously from a local 
Laplacian-type coupling (when a — * oo), to a global 
"mean- field" coupling (a = 0). Periodic boundary con- 
ditions and random initial conditions are assumed. 

A numerical diagnostic of CS is provided by the order 
parameter 



1 N 

z n = R n exp (2Tri(p n ) = — exp(27ria4 j) ) 



(2) 



where R n and ip n are the amplitude and angle, respec- 
tively, of a centroid phase vector, for a one-dimensional 
chain with periodic boundary conditions 10]. A time- 
averaged amplitude R is computed over an interval large 
enough to warrant that the asymptotic state has been 
achieved by the lattice. A CS state implies R = 1, for all 
the phase vectors gyrate together with the same phase. 
When the maps are completely non-synchronized, on the 
other hand, we get a pattern in which the site ampli- 
fy) 

tudes Xn are so spatially uncorrelated that they may be 
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FIG. 1: Synchronized regimes in the parameter plane of cou- 
pling strength e versus effective range a, for a lattice with 
TV = 21 coupled logistic maps. 



considered essentially as random variables, such that R 
vanishes. 

Figure ^ represents a portion of the coupling parame- 
ter plane (strength versus effective range) exhibiting re- 
gions for which the coupled maps have (i) synchronized 
chaotic orbits (CS); (ii) a transitional regime where CS is 
eventually attained after a transient characterized by in- 
termittent bursts of non-synchronized behavior; and (iii) 
completely non-synchronized orbits. These regions are 
bounded by the curves e*(a) and e c (a). The former cor- 
respond to the values of e* and a for which the aver- 
age order parameter magnitude R ceases to be equal to 
the unity. The latter curve was computed by means of 
the finite-time Lyapunov exponents, as will be explained 
later on in this communication. 

Let us fix the effective range parameter at an inter- 
mediate value, say a — 0.4, and decrease the coupling 
strength e from its maximum value, 1.0, to zero. For e 
large enough we have R = 1, or a chaotic CS state. When 
e = e c (0.4) w 0.69 this CS state starts to be interrupted 
by intermittent bursts of non-synchronized behavior, but 
eventually the stationary CS state is achieved again. The 
bursting becomes more frequent as the coupling strength 
is further decreased and, at e = e*(0.4) w 0.57, the order 
parameter vanishes and the lattice becomes completely 
non-synchronized, never to return to a CS state again. 
Hence, the interval e c < e < e* characterizes a transition 
region for which the intermittent bursting is a transient 
phenomenon. 

This scenario is robust and present for a wide portion of 
the coupling parameter plane. As we increase the value 
of the effective range parameter a from zero, it turns 
out that the interval characterizing a transition region 
is pushed towards higher values of the coupling strength 
e. We can understand this fact since the higher the ef- 
fective range a is, the closer we are to a locally coupled 
lattice, in which only the nearest neighbors contribute in 
a significant way, such that it becomes increasingly more 
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FIG. 2: Normalized histogram for the relative number of syn- 
chronization plateaus with a given length, for e = 1.0 and 
a = 1.147. 



difficult to have a stationary CS. We call a c ~ 0.940 the 
range parameter value for which the first critical curve 
reaches the upper limit given by e c (a c ) = 1. Accordingly, 
a* « 1.146 is the value where the second critical curve is 
such that e*(a*) = 1. Hence, for a > a*, we will not ob- 
serve a stationary CS state, irrespective of how strong the 
coupling may be, i.e. the intermittent bursting continues 
for an arbitrarily long time. 

Near the second critical curve in Fig. ^the average du- 
ration of the laminar regions obeys a power-law scaling 
< t >^ (a — a*) 7 , when a — > a* - , and with 7=1/2 
within the numerical accuracy. This suggests that the 
transition to synchronization occurs through a crisis |ll| . 
The rationale for this analogy identifies a burst between 
laminar regions as a chaotic transient resulting from the 
collision of a chaotic attractor with an unstable peri- 
odic orbit. The average chaotic transient length in one- 
dimensional maps, like the logistic map f(x) — ax(l — x) 
at a = 4, obeys an identical scaling and with the same 
critical exponent |l2j |. This crisis is mediated through 
an unstable-unstable pair bifurcation at those points at 
which there is such a collision [l3|. 

Another distinctive feature of the intermittent transi- 
tion to synchronization is the universal character of the 
statistical distribution of the laminar regions. In Figure 
121 we present a histogram for the number of the lami- 
nar regions with respect to their lengths. Two different 
regimes are highlighted in the insets: for small laminar 
regions the histogram is well fitted by a power-law r~ ro , 
with w k, 1.5; whereas the scaling is exponential e~ KT 
for large intervals, with k w 10 -3 . The presence of the 
3/2-scaling plus the exponential decay shows that the in- 
termittent transition to CS is of the on-off type ^4[ . We 
note that, in the neighborhood of the unstable-unstable 
pair bifurcation which mediates the crisis, a given map 
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FIG. 3: Normalized histogram of the largest transversal finite- 
time Lyapunov exponents for different values of the effective 
range a and e = 1.0. 



of the lattice chain receives random kicks through the 
coupling with its neighbors. These kicks, in the vicinity 
of the bifurcation point, are responsible for the observed 
on-off intcrmittcncy. In fact, the existence of two dis- 
tinct scalings in Fig. [21 with a shoulder indicates the 
presence of noise in the on-off intermittent scenario, with 
a crossover time related to the noise level 0] . 

In the transition region, trajectories off but very 
near the synchronization manifold experience intermit- 
tent bursting, such that numerical diagnostics using an 
insufficiently large time interval could erroneously point 
out CS 0| . On the other hand, this intermittent burst- 
ing is directly related to the shadowing breakdown of 
chaotic trajectories in the synchronization manifold of 
the coupled map lattice x\ = ■ ■ ■ = . The inter- 
mittent transition to CS initiates through an unstable- 
unstable pair bifurcation occurring on the synchroniza- 
tion manifold. The unstable periodic orbit undergoing 
this bifurcation, and all its infinite pre-images, increase 
their unstable dimension by one unit. Since there re- 
mains an infinitely large number of periodic orbits that 
did not undergo this increase, the synchronization man- 
ifold fails to be hyperbolic through unstable dimension 
variability 13]. 

This violates one of the mathematical requirements 
for hyperbolicity and shadowing, which is the continu- 
ous splitting of stable and unstable invariant subspaces 
along the attractor Q ■ If we start off but very near the 
synchronization manifold, after it has lost hyperbolicity 
through UDV, the corresponding trajectory approaches 
unstable orbits with different unstable dimensions. As a 
result, there are both time intervals for which the trajec- 
tory is, on average, approaching the manifold, and other 
intervals for which it is repelled, also on average. This be- 
havior is quantified by the largest finite-time transversal 
Lyapunov exponent, Aa(ri) 9]. If it is negative (posi- 
tive) a trajectory close to the CS state will be attracted 
(repelled) on average from it. The presence of UDV gen- 
erates random fluctuations of A2(n) around zero, in such 
a way that we can study its statistical distribution |18| . 
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FIG. 4: Average value and variance (inset) for the distribution 
of the finite-time Lyapunov exponents closest to zero versus 
a, when e = 1.0 



We get a numerical approximation for the probabil- 
ity distribution of this exponent, denoted as P(\2(n)) 1 
by considering a large number of trajectories of length 
n from initial conditions randomly chosen in the syn- 
chronization manifold. In Figure |31 we show some dis- 
tributions of time-50 exponents, obtained for different 
values of the effective range a, which here plays the role 
of a bifurcation parameter. The critical value for the 
onset of UDV is also the onset of intermittent transi- 
tion to CS, at a = a c , where an unstable- unstable pair 
bifurcation occurs. Numerically the onset of UDV was 
estimated through computing the value of a for which 
the fraction of positive finite-time Lyapunov exponents, 
<j>(n) = J °° P(\2(n))d\2, becomes non-zero, yielding the 
points on the curve e c (a) depicted in Figure ^ Hence, 
for a < a e , no shadowing breakdown via UDV is ex- 
pected, and the chaotic synchronized trajectories ob- 
tained through the numerical solution of Eq. (1) are 
expected to be shadowed over a longer time interval [(J, 
but which may be enough for practical purposes (e.g., 
when computing dimensions and entropies). 

The shape of the probability distributions in Fig. 
is Gaussian-like, with different variances a 2 , according 
to the value which a takes on. The Gaussian nature 
of P(\2(n)) is already expected on general grounds [Isj . 
and the distribution as a whole drifts toward positive val- 
ues of A2(n), as a increases. When a — a* the average 
< \2{n) > crosses zero and have a maximum UDV (Fig. 

since, at this point, < \2(n) > equals the infinite-time 
limit of A2(n), and UDV is more intense as the average 
of the time-n distribution crosses zero. This means that 
there are counterbalancing contributions from the nega- 
tive and positive finite-time exponents, as illustrated by 
the distribution for a — 1.14834^ a* shown in Fig. [3] 
The point at which < Aa(n) >= also marks the loss 
of transversal stability of the synchronization manifold, 
or a blowout bifurcation |l9j |. Accordingly, the points on 
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the curve e* (a) (Fig. ^) were also computed by imposing 
that < A2(n) >= 0, which furnished the same results as 
those obtained with help of the order parameter R, thus 
confirming the relationship between the loss of synchro- 
nization and the shadowing breakdown via UDV. 

In the vicinity of the point a = a* the shadowing times 
can be very short, and the validity of the computed nu- 
merical trajectories is doubtful. There results the time-n 
exponents (with n greater than the shadowing time), may 
suffer from similar shadowability problems, when tak- 
ing individually, as the chaotic trajectories themselves. 
However, in terms of the numerical diagnostics of UDV, 
we are actually interested in statistical properties of the 
time-n exponents, as their averages and variances. The 
former yields the point where UDV is the most intense, 
through a blowout instability, whereas the latter can be 
used to estimate shadowing times |5|. On the other hand, 
in some physically relevant cases, statistical quantities 
like these have been found to be meaningful despite un- 
stable dimension variability 20] . 

For a > a* , the relative number of positive time-n ex- 
ponents increases and we progressively return to a situa- 
tion where UDV is less pronounced, and better shadow- 
ing properties are expected. For example, when a = 2.5, 
Fig. ^indicates that almost all time-n exponents are pos- 
itive. This means that, even though the lattice trajecto- 
ries are far from being CS, they are nonetheless better 
shadowed than before. Figure0]also depicts the variance 
of the distributions, showing that it is roughly constant 
until a reaches the blowout value a*, after that the dis- 
tribution broadens up and the variance increases. 

In conclusion, we present numerical evidence that the 
intermittent transition to complete synchronization in 
a coupled map lattice is followed by a strong violation 
of hyperbolicity in the synchronization manifold. This 
breakdown begins with an unstable-unstable pair bifur- 
cation, through which an unstable periodic orbit and all 



its pre-images embedded in the synchronization manifold 
lose transversal stability. This is the typical way in which 
invariant manifolds lose transversal stability in riddling- 
blowout bifurcation sequences. As a consequence of hav- 
ing been left with an infinite number of transversely sta- 
ble orbits, the dimension of the unstable subspace varies 
from point to point along the synchronization manifold. 
There are severe consequences from the point of view of 
the validity of numerical trajectories obtained by com- 
puter simulation of high-dimensional extended systems, 
from which a coupled logistic map lattice has been chosen 
as a fundamental and representative example. The nu- 
merical trajectories may lose their validity, since they are 
no longer shadowed by true chaotic orbits for a reasonable 
time. The shadowing properties become worse as we ap- 
proach the point whereby the synchronization manifold 
loses transversal stability. The statistical properties of 
transverse finite-time Lyapunov exponents indicate that 
both the completely synchronized and completely unsyn- 
chronized maps have adequate shadowing properties, and 
the validity of numerical results within those realms is 
fairly well established. The shadowing problems occur 
just at the transitional regime, where laminar regions 
of chaotic synchronized behavior coexist with irregular 
bursts of non-synchronized activity. While our discus- 
sion was based on a specific coupled map lattice, the gen- 
eral features obtained here are rather typical for complex 
systems presenting regular and chaotic behavior in space 
and time. The connection between shadowing breakdown 
and intermittent transition to complete synchronization 
identified in this work is a common feature of complex 
systems. 
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